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Abstract 

Several physical applications of Lax equation require its general solution for 
generic Lax matrices and generic not necessarily diagonalizable initial conditions. 
In the present paper we complete the analysis started in [arXiv:0903.3771J on the 
integration of Lax equations with both generic Lax operators and generic initial 
conditions. We present a complete general integration formula holding true for 
any (diagonalizable or non diagonalizable) initial Lax matrix and give an original 
rigorous mathematical proof of its validity relying on no previously published results. 
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1 Introduction 



Lax equation appears in a variety of physical-mathematical problems which turn out to 
constitute integrable dynamical systems. 

An integration algorithm for generic matrix Lax equation was originally derived in 
the mathematical literature in [1], [2] and it was applied in the context of supergravity 
to cosmic billiards in [3], [1], [5] and to black-holes in [6], [7], [8] on the basis of their 3D 
description pioneered in [9]. 

The latter application to the case of black-holes revealed that the integration algorithm 
of [I], [2] did not cover the case of non-diagonalizable initial conditions. Such a situation 
occurs in particular when the initial Lax operator at time t = is nilpotent and this 
case is extremely relevant for Physics since it corresponds to extremal black- holes [9j. An 
extension of the integration algorithm to the case of nilpotent operators was presented 
in [8j. In the appendix of that paper it was also conjectured a general formula, which 
was verified for some nontrivial cases, that provides the integration of Lax equation for 
completely generic initial data. 

In the present paper we recall the general formula of [8] and we provide the rigorous 
mathematical proof that indeed it solves Lax equation with arbitrary initial conditions 
(diagonalizable or non diagonalizable) . Our proof is completely original and independent 
from the algorithm discussed in [I],[2]- 

2 The integration algorithm for generic Lax matrices 
and generic initial conditions 

Let us consider Lax equatioE0 



where L{t) denotes a time-dependent generic N xN matrix and remove any hypothesis on 
the nature of the initial conditions. We are interested in writing an integration algorithm 
for eq. (l2.ip which should hold true for generic initial matrices L{0) = Lq independently 
from the fact that Lq be diagonalizable or non diagonalizable, nilpotent or not, symmetric 
or non symmetric with respect to any definite or indefinite metric rj. 

It turns out that such an integration algorithm exists, is very simple and equally simple 
and elegant is the formal proof of its validity. 

Following [S] we first present the integration formula and then provide the mathemat- 
ical proof that it satisfies Lax equation. 

"'^ Hereafter, we denote by L> (i<) the upper (lower) triangular part including the diagonal of the 
matrix i, and L'^ is the transposed matrix. 
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L{t) + [L^it) - L^{t),Lit)] = 




1 



2.1 Integration formula 

The matrix elements Lpq{t) of the time-evolving Lax operator which coincides with the 
given Lq at time t = are constructed as follows [Sj: 



Lpq{t) 



1 



V 1 



Y,Y.'^r>k{t) {C{t)Lo),,mqe{t). (2.2) 



The building blocks appearing in eq. fl2.2p are defined as follows. We have 



\i+k 



Det 



V Ci,i(t) ... Ci,,-i{t) J 



l<k <i;2<i<N , 



(2.3) 



where the hats on the entries corresponding to the k-th row mean that such a row has 
been suppressed giving rise to a squared (i — 1) x (i — 1) matrix of which one can calculate 
the determinant. Similarly: 



f Ci,i(t) 



l)*+'=Det 



Ci,fc(t) 



\ Cj_l,l(t) . . . Cj_l_A;(t) . . . Ci-l^t) J 

1 < k <i ; 2 <i <l<i ; ^ii(t):=l, 



(2.4) 



where the hatted fc-th column is deleted just as in eq. fl2.3p it was deleted the k-th row. 
In the above formulae we have used the following definitions: 



-2tLo 



and 



Di(t) := Det 



Cit) 

\ Ci,i(t) ... Ci^t) ) 



2)o(t) :=1. 



(2.5) 



(2.6) 



2.2 The theorem and its proof 



In order to prove the integration formula (12. 2p we first recast it into another equivalent 
form which is not only the most convenient for the formal proof but also the simplest and 
best suited for computer implementation. 



2 



Theorem 2.1 The solution of Lax equation ^2. 1\) with generic initial condition L(0) = Lq 
is given by 

L(t) = Q{C)Lo (QiC))-' 



(2.7) 



where the N x N matrix Q(C(t)) is 



v/D.(t)2),_i(t) 

^ = Bet 

v/Di(t)Di-i(t) 



/ Ci,i(t) ... Ci,,_i(t) (Ci(t))i,, \ 



{21 



andC^t) = e"*^« (see, eq. /fO)). 

Proof 2.1 The matrix Q(C) satisfies three propertied 

X^{t) := Q(C)C2(t) = upper triangular , 
(A'<(t))~^ := Q{C) {C^t)y^ = lower triangular. 
{Xy{t)h = {X^{t))u 



(2.9) 
(2.10) 
(2.11) 



which are crucial in what fonows. Therefore, the matrix e = C^^t) admits the 
foUowing two different representations: 



-tLn 



-tLn 



X<{t) Q{C) 



(2.12) 
(2.13) 



resulting from eqs. (12. 9p and (]2.10p . Differentiating with respect to time t, one after the 
other equations f l2.12p . (12.131) and (12. 7p . then using them one can straightforwardly derive 
the following relations: 



s(c)|(S(c))- 



-1 



Q(C)| (Q(C))-^ = +L{t) + {X^{t)r' {jt^<(t)) 



dt 



m + 



Q{C) {Q{C)r\ L{t) 



(2.14) 
(2.15) 
(2.16) 



^The matrix {Xy{t))-- (|2.9p is upper triangular, since at j < i there is a coincidence of two columns 



in the determinant originating from 



m eq. 



The matrix (A'<(t))- -^ (|2.10p is lower triangular, 



since at j > i the all entries of the last column in the determinant originating from 
become equal to zero. 



in eq. ([^1^ 
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A simple inspection of eqs. (12 .141 - 12. 151) with the use of relation (12. lip leads to the conclusion 



that the quantity Q{C) | (Q(C)) is a traceless matrix which is expressed in terms of 
the Lax operator L{t) as follows 



Q(C)|(Q(C))-^ = L>(t)-L<(t) 



(2.17) 



then with this matrix Q{C) ^ {Q{C))~^ equation (12.161) reproduces Lax equation (12.11) 
and formulae (l2.7H2.8p give indeed its general solution for generic initial conditions Lq. 
This ends the proof of our proposition. (} 

Now, we present the useful expression of the inverse matrix (Q(C)) ^: 

1 



(Q(c))7/ 



( 



v/D^(t)D.-i(t) 



Det 



C.-i,i(t) 



(2.18) 



which follows from eq. (12.21) . It is a very simple exercise to verify that indeed 



(Q(C) (Q(C))-^) 



1 



VS)p(t) S)p_i(t)S),(t)D,_i(t) 



J] ^ $mpfc(t)Cfc,(t) aji,,(t) 



k=i e=i 



(2.19) 

since it is equal to zero if p 7^ g, due to a coincidence of two rows or columns in the 
determinants originating from (I2.8H2.18|) in eq. (12.191) . 

It is interesting to note one more property of the matrix Q{C) (12.81) : 



[2.20) 



which is obvious if one compares equations (12. 8p and (12.180 . When Lq is symmetric 
the same is true of C"^ and this implies that Q{t) is orthogonal, namely Q{t) G SO(N). 
Similarly when Lq is pseudo-orthogonal //-metric symmetric, Q{t) is pseudo-orthogonal, 
Q{t) G SO(p, N — P). Yet Q{t) exists in general also for non //-symmetric initial data. In 
this case Q(t) E qI(N). 



2.3 The Generalized Linear System 

Finally, based on the knowledge of the constructed general solution (I2.7H2.8I) with generic 
initial data Lq, we comment briefly on a way of a generalization (suggested by this solution 
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with the aim to reproduce it) of Kodama-Ye algorithm, which was based on the inverse 
scattering method apphed to the case of the diagonahzable Lq p]. 

A hnk towards the inverse scattering method is given by the following observation. 
Relations (12 .7^ and (12.171) corresponding to solution Q{C) (12.81) . being identically rewritten 
as follows: 



represent a generalized linear system for the matrix Q{C), whose consistency is provided 
by Lax equation (12. ip . This generalized linear system could be a starting point to apply 
the inverse scattering method to construct the solution for Q{C). But, there is a subtlety: 
up to now the inverse scattering method was always applied only when Lq is a diagonal 
or a diagonahzable matrix; it was never applied before to the case of non-diagonalizable 
Lq. Just the classical case of diagonal (diagonahzable) Lq was considered in [2] where the 
corresponding solution to the linear system was constructed. The very existence of our 
solution Q{C) (12. 8p for generic Lq gives an evidence that a generalization of the algorithm 
of [2j to the case of the generalized linear systems (12. 2 1112^221) with generic Lq should exist 
as well. It is indeed the case, and the generalization is rather straightforward. So it turns 
out that one can repeat all the steps of the algorithm, starting from the generalized linear 
system (I2.21f[2l22|) . in such a way that all the intermediate operations are formulated in 
purely matrix terms; thus they are insensitive to whether Lq is a diagonahzable or a non 
diagonahzable matrix: this process ends with the solution presented in (12.80 . 

3 Conclusions 

In the present paper we completed the analysis started in [8] on the integration of Lax 
equations with both generic Lax operators and generic initial conditions by presenting 
a simple, original rigorous proof of the validity of the integration algorithm proposed 
in [S]. The presented proof is constructive, and moreover it is interesting for its own 
sake since it opens a new deep insight into the structure of the corresponding integrable 
Toda-like systems which we plan to discuss elsewhere in connection with relevant physical 
problems like the classification and construction of Supergravity Black Hole solutions. 
We also clarified the relation between our new integration algorithm and the inverse 
scattering framework adopted by Kodama et al for the integration of Lax equation in the 
diagonahzable case. 
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L{t) Q{C) 



Q{C)Lq, 

- (L>(t) - L<(t)) Q{C), 





(2.22) 
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